Journal of Experimental Zoology Part B: Molecular and Developmental Evolution

Regenerative capacity varies widely across flatworms (Platyhelminthes). Whereas catenulids, microstomids and planarians can regenerate a complete head de novo, other flatworms cannot. This striking diversity raises a longstanding evolutionary question: does whole-body regeneration represent an ancestral trait that was subsequently lost in multiple lineages, or did it evolve convergently? Addressing this question requires comparative analyses of the molecular mechanisms underlying regeneration across phylogenetically diverse flatworms. Here, we focus on Wnt signaling, a deeply conserved regulator of antero-posterior (A-P) patterning and a central determinant of head-versus-tail identity during planarian regeneration, to establish a mechanistic framework for such comparisons. Although Wnt signaling has been studied extensively in planarians and parasitic neodermatans, its evolution and deployment in other flatworm clades remain poorly characterized. To address this gap, we characterized the complement of Wnt signaling components in two early-diverging flatworm clades, Catenulida and Macrostomorpha, with particular emphasis on expression and function in the catenulid Stenostomum brevipharyngium. Phylogenetic analyses reveal the ancient loss of six Wnt families and one secreted Frizzled-related protein (sFRP) family in the last common ancestor of flatworms, followed by additional lineage-specific gene losses and expansions. Moreover, several Wnt pathway components display markedly divergent expression patterns between catenulids and other flatworms, while functional analyses indicate corresponding differences in their regenerative deployment. Together, our findings reveal a dynamic evolutionary history of the flatworm Wnt signaling toolkit and establish a comparative framework for testing whether the molecular circuitry underlying head regeneration is ancestrally conserved or has evolved independently in distinct flatworm lineages.

Epithelial cells in Hydra perform an unusual combination of functions: they divide continuously like adult stem cells while simultaneously executing the complex physiological tasks of differentiated epithelia. This challenges the traditional distinction between proliferative stem cells and terminally differentiated tissue, raising the question of how a single cell type integrates these opposing roles. Using electron microscopy, we examined morphological characteristics that define the stem-like and differentiated states of Hydras ectodermal and endodermal epithelial cells. Stemness is reflected by nuclear characteristics of active proliferation, including extensive euchromatin, large nucleoli, and the presence of nuage. However, differentiated epithelial cells exhibit strong apical-basal polarity, various endomembrane compartments for endocytosis and transport, specialized secretion mechanisms, and basal muscle processes with dense-core vesicles implicated in hormonal communication. Cryofixation improved ultrastructure preservation, elucidating the pleiomorphic configurations of complex intracellular channel systems traditionally presenting as singular vacuoles. This may shed new light on possible functions of this compartment. Taken together, Hydra epithelial cells combine ancient stem cell traits with highly specialized differentiated functions. This multifunctionality provides insight into the cellular organization of early-branching animals and suggests that multifunctional epithelia may represent an ancestral condition preceding the strict segregation of stem and differentiated cell lineages in bilaterians.

Indomethacin and temperature shocks are both widely used to induce ephyra production in Cassiopea polyps, yet the combined effect of these stimuli on strobilation rate, ephyra viability, and polyp survivorship has not been systematically examined. We designed a full factorial experiment crossing five indomethacin concentrations (0, 10, 25, 50, and 75 {micro}M) with four temperature treatments (22, 25, 28, and 31{degrees}C) across three independent experimental trials, exposing a total of 288 individual Cassiopea xamachana polyps for 29 days. Overall, 207 of 288 polyps (71.9%) released at least one ephyra. Temperature was the stronger predictor of strobilation: higher temperatures substantially increased the probability of producing a healthy ephyra (log-odds: +0.327 per {degrees}C; {chi}{superscript 2} p = 6.5 x 10-5), while indomethacin had a smaller but significant positive effect (log-odds: +0.024 per {micro}M; {chi}{superscript 2} p = 0.035). The two stimuli acted independently without significant interaction on healthy ephyra production ({chi}{superscript 2} p = 0.87), although their interaction on total ephyra production (including unhealthy releases) was significant ({chi}{superscript 2} p = 1.0 x 10-9). Notably, low indomethacin concentrations (10 and 25 {micro}M) significantly increased the proportion of unhealthy ephyrae relative to seawater controls (paired t-tests: p = 0.035 and p = 0.006, respectively), while higher concentrations (50-75 {micro}M) did not. Polyp survivorship declined significantly at 28{degrees}C and 31{degrees}C but was unaffected by indomethacin. Trial replicate was the most statistically significant predictor of every outcome measured with the third trial producing near-universal strobilation but zero viable ephyrae, a result of an unidentified perturbation. These results indicate that temperature elevation to 28-31{degrees}C is the most reliable induction strategy for C. xamachana and that low indomethacin concentrations should be avoided when ephyra quality is paramount. We also note that cryptic batch-level variables can overwhelm controlled factors and should be addressed in future experimental designs.

Quantifying skeletal mineralization phenotypes in larval fish is complicated by the natural curvature of the notochord and by sample-to-sample variability in orientation, staining and imaging. Consequently, many studies rely on summary measures such as vertebral counts or total stain intensity. Here we present SCAMP (Spinal Calcification & Mineralization Profiler), an open-source, GUI-based Python tool that computationally straightens the curved notochord of Alizarin Red S-stained fish larvae and generates standardized mineralization profiles along the spinal axis. This approach reduces positional and shape variability, allowing direct, quantitative comparison of calcification patterns within and between experimental cohorts, without requiring programming expertise. We validate SCAMP using a zebrafish model of Pseudoxanthoma elasticum (abcc6aelu15/elu15), recovering genotype-specific differences in the intensity, extent and spatial distribution of ectopic calcification. Using SCAMP, we further show that inorganic pyrophosphate (PPi) supplementation of the medium suppresses ectopic notochord calcification, alters the anterior-posterior distribution of mineralized regions in homozygous mutants, and promotes mineralization at physiological vertebral sites. We also show that methylene blue, a routine antifungal additive in fish medium, reduces baseline calcification, with the most pronounced effects observed in heterozygous controls. SCAMP is freely available and has the potential to be adapted to other fish species used in skeletal and mineralization research.

Sea urchins are pivotal models in cellular and developmental biology, but their biphasic life cycle with an extended larval life limits the ability to study post-metamorphosis and adult characters. Here, we introduce a genomically-enabled model system, the tuxedo urchin Mespilia globulus, which has rapid access to late life stages in sea urchins. We describe how we cultured M. globulus in a landlocked marine facility, raised larvae under artificial conditions and closed their life cycle. We established the experimental tractability of M. globulus: we labelled transcripts by hybridization chain reaction (HCR), and knocked out pigmentation genes to produce albino larvae using CRISPR/Cas9. We generated chromosome-scale genome assemblies for two individuals representing both sexes and two color morphs (red and blue), and compared the organisation of the 21 chromosomes of M. globulus with that of other camarodont echinoid models. We determined that M. globulus showed a conservative gene repertoire lacking the gene family expansions seen in other camarodont sea urchins. We annotated the complement of genes associated with pigmentation, immune and nervous systems and profiled their expressions in tissues and organs. Finally, we surveyed sex-related regions in genomes using genome assemblies and resequencing data, finding no evidence of heteromorphic sex chromosomes in M. globulus. Our findings highlight the accessibility of this new sea urchin model for studying the metamorphosis and adult biology of sea urchins.

The nautilus (Nautilus sp.) is an early-branching cephalopods. It retains several conchiferan synapomorphies, including an external planispiral biomineralized shell. The shells complex structure allows it to withstand hydrostatic pressure, control buoyancy, and protect against external hazards. In this study, we comprehensively examined shell microstructures across different shell components and regions representing various ontogenetic stages in two adult museum shell specimens. We found that the nautilus shell is composed of five microstructural types (spherulitic, prismatic, nacreous, semi-prismatic, and irregularly oriented prismatic structures) organized into layered architectures within individual shell components and coordinated across the shell as an integrated system. Our observations highlight transitions between distinct microstructures within and across shell components and local variation within individual components such as the dorsal and ventral shell walls, suggesting that these patterns may contribute to shell strength and overall mechanical performance. Variation in caecum morphology suggests that this structure may be developmentally plastic and subject to relatively relaxed structural constraints. These findings show that the Nautilus shell is an integrated biomineral system in which diverse microstructures are organized across shell components to meet functional demands and provide the mechanical strength needed for survival.

Lactate dehydrogenase A (LDHA) is a key glycolytic enzyme that commonly exhibits altered expression in human diseases such as cancers and neurodegeneration, making it a valuable disease biomarker and putative therapeutic target. However, any treatment targeting LDHA will also disrupt normal metabolism, underscoring the need to investigate physiological consequences of inhibiting this enzyme. We previously established the fruit fly Drosophila melanogaster as a genetic model for studying LDH function in the context of growth, metabolism, and development. Here we expand upon those studies by investigating a serendipitous observation that Ldh mutant larvae exhibit diet-dependent lethality. Using a multiomic approach, we discovered this diet-dependent phenotype is independent of nutritional composition. Instead, Ldh mutant larvae are exercise intolerant and display reduced mobility, rendering mutant larvae sensitive to food consistency. Moreover, tissue-specific analysis reveals that LDH activity within muscle and peripheral glia are essential for larval viability raised on solid food. Intriguingly, these phenotypes mirror the pathophysiology of LDHA deficiency (Glycogen Storage Disease Type XI; GSD Type XI) in humans, where mild symptoms are exacerbated by physical exertion and environmental stress. Together, our findings further highlight the value of using Drosophila to explore the developmental and physiological consequences of Ldh inhibition.

Parental ageing can influence offspring through non-genetic mechanisms. The contribution of epigenetic inheritance parental effects still remains poorly understood. DNA methylation is a widespread regulator of gene expression that changes during development and ageing and can also act as a mediator of intergenerational effects. We tested whether age-related changes in parental DNA methylation are transmitted to offspring in the short-lived turquoise killifish (Nothobranchius furzeri, Cyprinodontiformes). Using reduced-representation sequencing, we quantified genome-wide DNA methylation and examined methylation dynamics at individual loci. The overall proportion of methylated CpG sites increased during early ageing but declined at later ages, revealing a non-linear trajectory with substantial among-individual variation. Despite these age-related changes, we found no evidence that parental methylation patterns were transmitted to offspring, either at the genome-wide level or at individual loci. Our findings indicate that although DNA methylation undergoes pronounced age-dependent remodelling in adult killifish, these changes are not detectably inherited by the next generation. We discuss these results in the context of epigenetic inheritance and ageing in vertebrate model systems.

Although corals are highly regenerative, some colonies in reef ecosystems completely recover from sublethal damage while other colonies exhibit partial mortality to similar injuries. To understand factors that might naturally curtail regenerative ability, we experimentally wounded small colonies in three coral genera (Acropora, Pocillopora, Porites) by mimicking natural corallivory using scraping (tissue and skeletal damage) or airbrushing (deep tissue removal with no skeletal injury). We found all scraped wounds regenerated rapidly in Acropora and Porites, while Pocillopora fragments frequently retained open lesions. In stark contrast, airbrushing resulted in algal colonization and delayed tissue healing and regeneration across all corals. Detailed cellular analysis of Porites wounds revealed two general phases comprising tissue regeneration: a healing phase defined by rapid coverage of bare skeleton with coenosarc, pigment cells and gastrodermal reformation, and then a second phase lasting one week ending in polyp regeneration. Red fluorescence appeared transiently in scrape wounds but persisted in tissue at the wound margin surrounding algae in airbrush wounds, suggesting that algal occupation of the wound bed inhibits coenosarc healing. Lastly, histological cross-sections of healing airbrush wounds in Porites revealed progressive loss of deep tissue leading to skeletal breakdown beneath the wound. Together, our results demonstrate the biphasic nature of tissue regeneration in colonial corals and provide a framework for understanding how biotic factors impact tissue repair and regeneration in nature.

Dauer formation and L1 arrest are stress-responsive developmental strategies that enable Caenorhabditis elegans to survive unfavorable conditions. These responses are regulated by environmental cues, including food availability and pheromone signals that communicate population density. However, how dauer entry, L1 arrest, and density-dependent signaling collectively influence long-term population dynamics remains poorly understood. In this study, populations of daf-22 and daf-16 mutants with impaired dauer formation, starvation arrest, and pheromone signaling were compared under control and starvation stress. Measurements of developmental stages were used to evaluate how genotype influenced population growth, developmental stage composition, starvation response, and recovery over time. This population-level approach links individual developmental decisions and inter-organismal communication to broader patterns of persistence and population change. The results show that daf-16 and daf-22 mutant populations differed from wild type in their recovery ability following nutrient deprivation as well as in the stage distributions within the population. These findings suggest that dauer signaling contributes broadly to population persistence by coordinating developmental arrest, reproduction, survival, and recovery. Overall, this work supports the interpretation of dauer formation as a larger population level survival program rather than a single isolated developmental outcome.

Drosophila imaginal discs are the favoured system to study organ size determination. The intrinsic mechanism of organ size determination remains elusive despite decades of research into the cellular nature, mechanics and genetics of growth in the wing imaginal discs. Here we propose a new model whose basis is growth rate dependent-negative feedback without size being sensed. We show that alterations in size before the growth period are not corrected for, but alterations in growth rate during the start of the growth period induce size-restoring feedback.

The influences of the gut microbiome on animal physiology are well-documented, yet the developmental timing of microbial colonization and its long-term consequences remain poorly understood. In this study, we investigated how the timing of bacterial colonization during development affects transcriptional programming and phenotypic outcomes in adult Drosophila melanogaster reared on a common, rich diet. Using RNA-seq analysis on whole flies colonized either as newly hatched larvae or as newly eclosed adults, we observed minor but distinct transcriptional responses dependent on when flies were colonized. Both embryonic and adult colonization were associated with ~ 25 to ~ 200 differentially expressed genes compared to axenic controls, with the majority upregulated and enriched for immune-response genes, suggesting that colonization establishes a broader immune competence. However, only 10 genes showed persistent differential expression that was not normalized by introducing bacteria to adult axenic flies, including mitochondrial genes, the adipokinetic hormone (Adh), and a putative secreted neuropeptide. Overall, these findings suggest that Drosophila development on a rich diet is largely robust to the timing of bacterial colonization but that certain metabolic effects may occur.

The physiology of ectotherms can change substantially during acclimation to changing environmental temperature. The role of transcription in acclimation responses has been well-established, but it remains unclear whether transcriptional regulation generally reflects abrupt changes after surpassing temperature thresholds, or whether transcript abundance is a relatively monotonic, continuous function of acclimation temperature. In this study we exposed adult male Drosophila subobscura flies to four different 96-hour acclimation treatments at temperatures that were not acutely stressful but ranged from relatively cold (10{degrees}C) to relatively warm (27{degrees}C) with respect to standard rearing conditions. Transcriptome sequencing of whole-body homogenates (mRNAseq) revealed a massive, transcriptome-wide response across acclimation temperatures, with a marked overrepresentation of genes that were continuously and monotonically up- and down-regulated in response to increasing acclimation temperature. Though some genes showed more complex relationships consistent with putative threshold responses, a high percentage of the differentially expressed transcriptome (42%) showed continuous and strictly monotonic relationships. Functional enrichment suggested continuous up-regulation of spermatogenesis-related transcripts with increasing temperature and continuous up-regulation of oxidative phosphorylation-related transcripts with decreasing temperature, illustrating contrasting patterns consistent with previous studies of thermal sensitivity of male reproduction and metabolic compensation in the cold. Thus, continuous thermal sensitivity of transcription is a hallmark of acclimation in D. subobscura that likely underlies the continuous thermal sensitivity of downstream physiological processes. We also provide evidence for shared transcriptomic responses across short-term acclimation (this study) vs. published results for long-term, developmental acclimation.

In a wide variety of animals, developmental crowding results in adults with smaller bodies. The crowding effect on body size in Drosophila melanogaster is canonically attributed to heightened competition for nutrition. However, whether other consequences of crowding also contribute to its effect on size remains an open question. We tested the relative contributions of nutritional competition, oxygen availability, and larval-generated metabolites to the crowding effect on size. We found that while nutrition explains most of the variation in body size due to crowding, oxygen also contributes in a sex- and nutrition-dependent manner. We found no evidence that larval-generated chemicals affect body size. These data confirm a widely suspected but untested role of nutrition in producing the crowding effect on size in D. melanogaster, while revealing an unexpected role of oxygen, and raise the possibility that behavior may be a mediator of density-dependent plasticity. Research HighlightsWe found that both nutrition and oxygen mediate the crowding effect on size in Drosophila melanogaster.

The sea star Patiria miniata is a widely used and powerful model organism for cell, developmental, and reproductive biology, but it has lacked genetic tools for expressing transgenes or endogenously tagging proteins. We developed a protocol to endogenously tag a broadly expressed actin gene and to express additional fluorescent markers from the same locus, using CRISPR/Cas9 genome editing. We also identified and isolated a promoter sequence of this actin gene which drives expression of transgenes. This promoter and transgene cassette can be introduced via a plasmid into the genome and persist through metamorphosis into the juvenile stage. Robust methods to induce metamorphosis that result in healthy juveniles are essential for developing stable transgenic lines and have been lacking in the field. Here we report a fast and efficient approach to induce the metamorphosis of larvae into healthy juveniles by introducing surf clam shells. Thus, we present a reliable method to generate both CRISPR/Cas9 knock-in and plasmid-integrated transgenic juvenile P. miniata, enabling future research on their fascinating biology, including regeneration, oogonial stem cells, metamorphosis, and more.

Host-associated microbiomes play vital roles in organismal health, ecological interactions, and evolution, yet little is known about how microbial communities shift during the transition from aquatic to increasingly terrestrial habitats. Here, we characterize the skin microbiomes of three combtooth blenny species (Blenniella paula, Praealticus labrovittatus, and Alticus arnoldorum) that occupy distinct positions along the intertidal gradient in Guam--from fully subtidal (and exclusively aquatic) to intertidal (amphibious) and supratidal environments (exclusively terrestrial). Using 16S rRNA amplicon sequencing, we compared skin-associated bacterial communities with those in surrounding seawater and substrate biofilms to assess habitat influences on microbiome structure. Skin microbiomes were distinct from environmental microbial communities, indicating strong ecological filtering by the host. The divergence between skin and substrate microbiomes in the three species parallels their distribution along progressively higher zones of the intertidal gradient. The most divergent skin microbiome was that of the supratidal fish A. arnoldorum, characterized by higher Gammaproteobacteria abundance and enrichment of epiphytic and mucus-associated taxa. Across all species, we identified 32 microbial orders significantly enriched on the skin relative to environmental samples, including taxa commonly associated with fish mucosa (e.g., Vibrio, Alteromonas, Cetobacterium) and others rarely reported in aquatic marine fish (e.g., Rubritalea, Granulosicoccus). Several rare taxa with potential pathogenicity were also detected at low abundances. Together, these findings suggest that habitat-specific selective pressures strongly shape fish skin microbiomes along subtidal (aquatic) to supratidal (terrestrial) habitats and suggest that microbial symbionts may contribute to the ecological and physiological adaptations enabling amphibious lifestyles. This study provides the first comparative assessment of skin microbiome divergence across amphibious fish species along an intertidal gradient and offers a framework for predicting microbiome responses to environmental change.

Multiple hypotheses have been formulated to explain differences in tissue repair ability across vertebrates. One hypothesis posits that the accessibility of actively cycling stromal cells within uninjured tissue confers access to a proliferative population in response to tissue damage. This hypothesis further suggests that animals with an indeterminate growth mode possess an actively cycling cell population necessary for growth that can be readily accessed for tissue regeneration. Moreover, the absence of an actively cycling population in connective tissue provides a mechanism that restricts regeneration in animals with determinate growth whose cells are refractory to cell cycle progression and proliferation to produce new tissue for morphogenesis. Here, we explore this paradigm using an EdU-BrdU pulse chase strategy in four different vertebrate species: two with determinate (Acomys dimidiatus and Mus musculus) and two with indeterminate modes of growth (Danio rerio and Ambystoma mexicanum). We find that although indeterminate growers do possess a small population of actively cycling cells, this population does not contribute to regeneration. Moreover, we found that while Acomys does not possess a population of actively cycling stromal cells, cells re-enter the cell cycle de novo in these animals to contribute to regeneration. Furthermore, testing this hypothesis allowed us to ask whether tissue injury could stimulate cell cycle re-entry - a so-called primed state - in cells at distance from the injury site in these four species and we did not find evidence of such priming in stromal or epidermal tissue. HighlightsO_LICell cycle re-entry is a common response to injury in regenerative and non-regenerative vertebrates that is independent of actively cycling stromal cells in uninjured connective tissue C_LIO_LIActively cycling cells do not contribute to regenerative healing in spiny mice, axolotls or zebrafish. C_LIO_LIOur data do not support systemic cell cycle activation in response to injury. C_LI

RNA interference (RNAi) is a widely exploited reverse-genetics tool with potential uses for disease control. Successful RNAi has been reported in trematode-vectoring snails, but the composition of RNAi effector-encoding gene complements, a key driver for RNAi efficiency, remain unstudied in these species. Using bioinformatics and comparative genomics, we searched for orthologues of 115 RNAi effector sequences in genomes or transcriptomes of four snail vectors: Biomphalaria glabrata, B. pfeifferi, Bulinus truncatus, and Lymnaea staginalis. Gene expression patterns of selected RNAi effectors were then examined across developmental stages and tissues of the model B. glabrata snail. At least 74 RNAi-related proteins were conserved across all four species, including core components known to be essential for gene silencing. Classical systemic RNAi-deficient (SID) genes that facilitate systemic RNAi in other systems were absent, suggesting that alternative pathways may compensate for dsRNA uptake and transport. Core effectors of secondary RNAi amplification and heritable RNAi were not detected. Expressions of Dicer-1, Argonaute-2, and the exonuclease Eri-1 did not vary significantly with snail size. A putative RNAi-inhibiting Staufen orthologue showed elevated expression in the ovotestis, while another putative cholesterol-interacting gene was overexpressed in the trunk tissue and may partly contribute to RNAi import. Altogether, our results present the most comprehensive overview of RNAi pathway effectors in major intermediate snail hosts for trematodes. The findings underscore the likely broad potential for RNAi use in trematode intermediate hosts as an experimental tool and potential control method.

Exposure to the potassium channel blocker barium chloride (BaCl2) causes head degeneration in Dugesia japonica flatworms, followed by regeneration of BaCl2-insensitive heads, offering a unique model for studying transcriptional resilience to novel stress. We performed RNA sequencing on individual planaria to investigate different transcriptional solutions to the BaCl2 challenge, and how regeneration history and social environment shape transcriptomic responses to BaCl2. We identified a robust transcriptional strategy and a potential sub-strategy for enabling BaCl2-insensitive head formation. Moreover, we observed pronounced transcriptional differences between untreated worms regenerating from tail fission fragments (tail-regenerated), and untreated full-sized worms that did not fission during the experiment (intact controls), highlighting the lasting impact of regeneration history. Relative to controls, tail-regenerated worms upregulated neurodevelopmental and morphogenetic programs, while downregulating mitochondrial transport and stress-response pathways. Relative to intact controls, BaCl2-exposed regenerates upregulated ion transport, metabolic, cell cycle, and inflammatory pathways, while downregulating neuronal signaling, ion homeostasis, morphogenesis, and tissue repair programs. Comparison of BaCl2-exposed isolated and BaCl2-exposed group-housed worms revealed minimal transcriptional divergence between social conditions. These findings underscore the complex interplay between regeneration, chemical stress, and social context in shaping gene expression.

Long-term persistence of managed species will depend, in part, on whether the species harbors the physiological or genetic potential to adjust to warming temperatures, and whether relevant genetic variation is modified by management practices. The critically endangered Delta Smelt (Hypomesus transpacificus) is intensively managed, but little is known about the presence of genetic variation for resistance to elevated temperature, which will be important to maintain for their persistence in a rapidly warming future. Using a pedigree and whole genome sequencing data, we characterized the genetic variation and genomic architecture for CTMax (as a metric of upper thermal tolerance) across control and elevated rearing temperatures, alongside covarying traits (body size, degree of hatchery ancestry). Warmer rearing temperatures increased CTMax through acclimation but also resulted in reduced additive genetic variation for the trait, which could constrain adaptation under thermal stress. We found that larger fish had reduced CTMax, although this effect was diminished at elevated temperatures. We observed modest heritability for CTMax at rearing temperatures of 15{degrees}C and 18{degrees}C (0.26 and 0.16, respectively), but only a limited number of loci were identified that had consistent effects on CTMax across rearing temperatures. Instead, the genomic basis of thermal tolerance was highly dependent on rearing temperature (many loci detected with a GxE effect). The influence of domestication selection was indicated by changes in allele frequency, and divergence in upper thermal tolerance and plasticity, between low and high hatchery ancestry groups. Minimal overlap between loci associated with domestication and CTMax suggests that these traits possess separate genetic underpinnings. Knowledge of genetic variation supporting ecologically relevant physiological variation may be useful for refuge management and may inform supplementation in an ever-warming environment.